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Structure-induced features of the wave functions for the quantum systems on complex networks 
are discussed in this paper. For a quantum system on a network, the state corresponding to the 
eigenvalue close to the center of the spectrum is used as the representative state to display the 
impacts of the structure on the wave functions. We consider the Erdos-Renyi, the WS small world 
and the growing randomly network (GRN) models. It is found that the probability distribution 
functions (PDF) of the representative state's components can be described with a power-law with an 
exponential cutoff in a unified way. For Erdos-Renyi networks, with the increase of the connectivity 
probability per the PDF turns from power-law-dominated to exponential-dominated functions. For 
the WS networks in a special region of the rewiring probability pr £ (0,0.2), where this model can 
capture the features of real world networks, and the GRN networks, the PDFs obey almost a perfect 
power-law. These characteristics can be used as the structure measurements of complex networks. 
They can also provide useful information on dynamical processes on complex networks. 
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The structures of complex networks can induce non- 
trivial features to the physical processes occurring on 
them. Typical topics include the epidemic spreading on 
networks P) S 13 > the synchronization of coupling oscil- 
lators on networks 0, IS 0, Q , the response of complex 
networks to stimuli [Sj, |9( and so on. In this paper, we 
consider the impacts of the network structures on the 
wave functions of quantum systems. 

Anderson transition tells us that disorder structures 
can induce a transition from extended to localized states 
[lO', "ll] . Quantum systems with quasi-period structures 
will be in an intermediate state, which can be described 
with critical wave functions 0, ^3 ■ The wave function 
for a localized state decreases exponentially with the dis- 
tance from its center, while the critical wave function 
obeys a power-law with respect to the distance. Com- 
plex networks have nontrivial structures rather than reg- 
ular and complete disorder ones. In literature the spec- 
tra density function and the time series analysis methods 
are used to reveal structure-induced features 
from the spectra HI E BB IH B 113 and the cor- 
responding eigenvectors |24ll25l| of complex networks. In 
this paper we try to capture the special characteristics of 
complex networks from the wave functions of the quan- 
tum systems on them. These wave functions can provide 
us useful information on the corresponding classical sys- 
tems, i.e., the probability distribution function (PDF) of 
packets walking on the networks at steady states [2(^ . 
Consequently, the PDF can shed lig ht on the dynam- 
ical processes as the traffic flow |23, |2^, the epidemic 
spreadiiig Q, |^ I3 and the synchronization on networks 

mm 

We consider an undirected complex network of N cou- 



pling identical oscillators. Denote the adjacent matrix of 
this network with A, whose element Aij is 1 and if the 
nodes i and j are connected and disconnected, respec- 
tively. The Hamiltonian of this quantum system reads. 



N ^ N 

H ho{Xn, Pn) + 2^ 

n—1 m^n 



V{x„i,Xn), (1) 



where ho{xn,Pn) is the Hamiltonian of the oscillator n 
and V{xm,Xn) the coupling potential between the oscil- 
lators m and n. Denote the site energy and the corre- 
sponding state of each oscillator with Bq and \ ipa), respec- 
tively. The elements of iJread, 



^ (^ipo{Xm) ho{Xm,Pm)\^o{Xn) 
+ A,nn ■ {iPa(xm) \ViXm,Xn) |<Po(x„) ) 



(2) 
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Assigning Eq = and Vmn = 1, we have H — A. Hence, 
the special features of the wave functions of the consid- 
ered quantum system can be obtained from the eigenvec- 
tors of the adjacent matrix A. 

The quantum systems considered in literature have 
deterministic structures in real world Euclidean space, 
which lead to nature definitions of the concepts as the lo- 
calized, intermediate and extended states of the quantum 
systems. These concepts are also extended to capture 
the features of quantum systems on small-world networks 
based on regular lattices, where the long-range edges can 
be regarded in a certain degree as perturbations to the 
regular lattices Obviously, these concepts are in- 

valid for general complex networks without deterministic 
structures in Euclidean space. Herein, we consider the 
probability distribution functions (PDF) of the values of 
the components in the eigenvectors. 
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The eigenvector corresponding to the special eigen- 
value close to the center of the spectrum for a net- 
work, denoted with E,,. is employed as the represen- 
tative state to illustrate the features of the consid- 
ered quantum system. Denote the representative state 
with r = (ri, r2, • • • , Tat), where is the mth compo- 
nent. Divide the range the probability values of p = 

|iV • |rm|^ , m = 1, 2, 3, • • • , distribute into M bins. 

Reckon the numbers of the values falling in the bins, re- 
spectively. The PDF of the values of the components can 
be approximated as. 



(3) 



where Kt is the number of the values falling in the tth 
bin. The size of the bins can be chosen as a fraction of the 



variance, i.e., Sb, 



N 



In the calculations 



the parameter J is assigned 10. It is found that the PDF 
for different networks can be described with a power-law 
with an exponential cutoff, which reads. 



P(r) oc T-"cxp(-/3T), 

T = t- (pmax - Pmin)/M. 



(4) 



For a = and /? = 0, the PDF will degenerate to expo- 
nential and power-law functions, respectively. 
The probability moment (PM) defined as. 



(5) 



m=l 



is also used to measure the global extent of localization. 
For a perfect extended state we have Mq{Ec) = f^l-i , 
while for a state strongly localized on one node the PM 
tends to 1. Generally, PM should be in the range of 

As an extreme condition, wc consider networks con- 
structed with the Erdos-Renyi (ER) network model. An 
ER network can be obtained just by connecting each 
pair of N nodes with the connectivity probability per- 
There is a critical point of the connectivity probability 
Per, denoted with pc, when per > Pc almost all the 
nodes tend to form a giant cluster. As shown in figure 1, 
with the increase of per- the value of the parameter /3 
turns from negative to positive. At the transition point 
Per = jf>Pc, we have /3 = 0.007 » 0. That is, the PDF 
is a perfect power-law function. When the connectivity 
probability Per = jf, P increases rapidly to 0.31 while a 
decreases to 0.87. The PDF is exponential-dominated. 

In the WS small-world model the nontrivial features 
are introduced into the regular lattices just by rewiring 
with a certain probability Pr the end of each edge to an- 
other randomly selected node. In the rewiring procedure 
double edges and self-edges are forbidden. In the regu- 
lar lattices each node is connected with d right-handed 
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FIG. 1: Typical P{T) results for the quantum systems on 
the ER networks. The circles and solid lines are the PDF 
and the fitting results, respectively. The power-law with an 
exponential cutoff presented in Eq.(4) can capture the char- 
acteristics of P(T) very well. The value of the parameter (5 
increases from negative to positive. At the transition point 
Per = > Pc, we have (3 = 0.007 « 0, that is, the PDF is a 
perfect power-law function. When the connectivity probabil- 
ity psj? = /3 increases rapidly to 0.31 while a decreases to 
0.87. Its PDF tends to be exponential-dominated. 
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FIG. 2: TypicalP(r) results for the quantum systems on the 
WS networks. The parameters {N, d) = (3000, 2). The circles 
and solid lines are the PDF and the fitting results, respec- 
tively. The power-law with an exponential cutoff presented in 
Eq.(4) can capture the characteristics of P{T) very well. 



nodes. Figure 2 presents several typical P(T) results for 
the quantum systems on the WS networks with differ- 
ent values of rewiring probability. From figure 3 we can 
find that in the special range of pr S (0, 0.2), where the 
model can capture the features of real world networks, 
the values of /3 are in the range of [—0.06,0.1], i.e., the 
PDFs obey almost a perfect power-law. The correspond- 
ing values of a are significant larger than that in the other 
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FIG. 3: (Color online) The values of {a, (3) for the quantum 
systems on the constructed WS networks. For each rewiring 
probability pr ten simulated results are presented with open 
circles, whose averages are shown with the red solid lines. 
{N,d) = (3000,2). In the interested region of Pr e (0,0.2) 
the value of (3 is in the range of [—0.06, 0.1]. The PDFs obey 
almost a perfect power-law. The corresponding values of a 
are much larger than that in the other regions. 



FIG. 5: Typical P(r)results for quantum systems on the 
regular lattices. The power-law with an exponential cut- 
off presented in Eq.(4) can capture the characteristics of 
P(T) very well. The values of (a, /3) decrease rapidly to 
(0.96, —1.307). The quantum systems on these regular lat- 
tices are significantly delocalized compared with that on WS 
networks (the values of the probability p distribute homoge- 
nously in a much narrower interval). The circles and the 
solid/dashed lines are the PDFs and the fitting results, re- 
spectively. (iV,d) = (3000,2). 
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FIG. 4: (Color online) The PM for the quantum systems on 
the WS networks with different values of rewiring probability. 
For each rewiring probability pr ten simulated results are pre- 
sented with open circles, whose averages are shown with the 
red solid lines. (iV, d) = (3000, 2). In the region pr € (0, 0.2), 
the PM values are much larger than that in the other regions. 
The quantum states on these networks are localized much 
more than that on the other networks. 



regions. Figure 4 presents explicitly the higher extent of 
global localization of the states for the quantum systems 
in the region pr £ (0, 0.2). For WS networks with d = 5, 
we can obtain similar results. 

For the conditioner — 0, the constructed WS networks 
degenerate to one-dimensional regular lattices. The val- 
ues of (a,/?) decrease rapidly to (0.96,-1.307). The 



FIG. 6; TypicalP(r) results for the quantum systems on the 
CRN networks. (iV,m) = (3000,1). The circles and solid 
lines are the PDF and the fitting results, respectively. The 
power- law with an exponential cutoff presented in Eq.(4) can 
capture the characteristics of P{T) very well. 



quantum systems on these regular lattices are signifi- 
cantly delocalized compared with that on WS networks 
(the values of the probability p distribute homogenously 
in a much narrower interval). 

The growing randomly network (GRN) model simu- 
lates a kind of preferential attachment for new added 
nodes during the growing of the networks. Starting 
from several connected nodes as a seed, at time t -t- 1 
a new node is linked to m existing nodes at time t, 
denoted as {sj|i=l,2,---,m}, with the probability 
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FIG. 7: (Color online) The values of (q, /3) for quantum sys- 
tems on the constructed GRN networks. The values of (3 are 
basically in the region of [-0.05,0.05]. The PDFs for the 
quantum systems on these GRN networks obey almost a per- 
fect power-law. {N, m) = (3000, 1) 
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FIG. 8: (Color online) The PM for the quantum systems on 
the constructed GRN networks. Comparison with the PM 
results for WS small-world networks shows that the quantum 
systems on GRN networks are much more localized than that 
on WS small- world networks. {N,m) = (3000, 1). 



p{s\) oc k(s\Y,Q < 6 < 1. k(s\) is the degree of the 
node s\. = \ leads to the BA scale-free model. Several 
typical results are shown in Figure 6. The values of (a, (3) 
for quantum systems on the constructed GRN networks 
are shown in figure 7. The values of /3 are basically in the 
region of [—0.05,0.05]. Consequently, the PDFs for the 
quantum systems on these GRN networks obey almost a 
perfect power-law. The global extents of localization are 
displayed in figure 8. 

Figure 9 presents the relation of a versus (3. In this 
scheme each point corresponds to a network. In the re- 
gion of /3 « the PDF obeys perfect power-law, while a 



FIG. 9: (Color online) The relation of a versus /3. For same 
values of 13, the a values for quantum systems on the GRN 
networks are generally smaller than that on the WS networks. 
There are much more nodes having large values of p in GRN 
networks. The ER networks with per less than or near the 
critical value Pc {per < jj) fall into the region of GRN net- 
works, while that with Per > jf falls in the region of WS 
networks. The network with per ~ jj falls in the extension 
of the WS network region. Its large value of /3 = 0.31 tells us 
that its PDF is exponential-dominated. 

significant deviation of /3 from reveals the exponential- 
dominated behavior of the PDF. For same values of /3, 
the corresponding a values for the quantum systems on 
the GRN networks are generally smaller than that on the 
WS networks. There are much more nodes having large 
values of p in GRN networks. This higher extent of local- 
ization can be found by comparison the PM values. The 
ER networks with per less than or near the critical value 
Pc [per. < j^) fall into the region of GRN networks, while 
that with Per > jf falls in the region of WS networks. 
The network with per — jf falls in the extension of the 
WS network region. Its large value of f] — 0.31 tells us 
that its PDF is exponential-dominated. 

In summary, the states of the quantum systems on 
the WS, GRN and ER networks obey a unified power- 
law with an exponential cutoff. The WS networks with 
Pr € (0,0.2), where the WS model can capture the fea- 
tures of real world networks, and the GRN networks obey 
almost a perfect power-law. The values of the power-law 
exponent a for GRN networks are much smaller than 
that for the WS networks. With the increase of per the 
PDF for ER networks tends from almost perfect power- 
law to exponential-dominated functions. These findings 
are consistent with the global extent of localizations for 
the considered networks, which can be displayed explic- 
itly with the PM values. 

The two parameters introduced in this paper, denoted 
with a and /3, are determined by the structures of the 
considered networks. Consequently, they can be em- 
ployed as a measurement of the structure characteristics. 
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The PDF can be used to determine the distribution of 
traffic packets on networks at steady states. This infor- 
mation may shed fight on the dynamical processes on 
networks, such as the traffic flow, the traffic congestion 
and the epidemic spreading processes on complex net- 
works. 
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